Electron transport driven by a chemical potential difference 
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Based on Bhatnagar-Gross-Krook equation coupled with Maxwell equation, we investigate 
the spatial dependence of a chemical, an electrostatic and an electrochmeical potentials inside 
a specimen connected with reservoirs. We also confirm that a gap of the chemical potential 
between at a connection point is negligible. 
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Study of the electron transport has a long history^ 
and is one of the typical problems in the nonequilibrium 
statistical physics. In these days there has been an in- 
creasing interest in understanding the electron transport 
by the successful fabricated devices.^ In the nonequil- 
brium transport theory, it is noted that electron currents 
and heat currents are related with the temperature gra- 
dient, the chemical potential gradient and the electric 
field. For these thermoelectric relations, Wiedemann- 
Franz law^ and Seebeck effect'' are well known. On the 
other hand, it is interesting to investigate how the local 
thermodynamical quantities are distributed in the con- 
ductance. In this leter, we focus on the spatial depen- 
dence of the thermodynamical quantities in the steady 
electric current which is driven by the chemical potential 
difference. 

Let us consider the electronic device which is com- 
posed of a specimen connected through reservoirs. The 
electric current is driven by the difference of the chemi- 
cal potential in both sides of the reservoirs. We assume 
that the system is at the isothermal low temperature and 
any effects of the heat current do not exist. Furthermore, 
Joule heat caused by the electric current assumed to be 
negligible. At the connection point between the speci- 
men and the reservoir, we need to be careful about the 
discontinuity of the thermodynamical quantities.^' ^ The 
schematic picture of the device is given by Fig. 1. The 
length of the specimen is larger than the inelastic scat- 
tering length and its specimen is directly connected with 
two equilibrium reservoirs. Here, we do not consider the 
collimator effect at the connection point. Each reser- 
voir has the chemical potentials /i^ and which keep 
constant values during the steady electric current. 

Although the interactions between the electrons play 
an important role for strongly correlated electron sys- 
tems, we can use the nearly free electron model to ex- 
plain the electric current in the metal. In this letter, 
we treat the electrons as semi-classical particles which 
are represented by a nonequilibrium distribution func- 
tion. To describe such electron transports, we adopt BGK 
(Bhatnagar-Gross-Krook) equation^ which is Boltzmann 
equation with a relaxation time approximation. Further- 
more, we couple BGK equation with classical Maxwell 
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equation. Under this setup, we investigate the spatial de- 
pendence of the chemical, electrostaic and electrochem- 
ical potentials in the specimen, and the discontinuity of 
the chemical potential at the connection point. 

We start with BGK equation to treat the electron 
transport. Here, we are interested in the x-directional 
spatial dependence of the thermodynamical quantities, 
so that we consider the averaged behavior of the elec- 
tric current in the cross section. The steady-state BGK 
equation for the distribution function f{p) is 



(1) 



where — e is the charge of the electron, E is the x- 
component of the electric field, Px is the a;-component 
of the linear momentum and the collison term /[/] is 
given by 



m = 
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where r is the relaxation time and — 1/ (e'^^'^P^^^^ + 1) 
is Fermi distribution function with the equilibrium chem- 
ical potential /x, and e(p) — /2m is the kinetic en- 
ergy of the electron. The electric field has a contribution 
from the chemical potential difference and an internal 
contribution from the induced charge density fluctua- 
tions. Here we do not consider the magnetic field induced 
by electric current. When the distribution function / is 
slightly different from the equilibrium Fermi distribution 
function,"^ we obtain the solution 



f = f-f{l-f)ePrVx[E + 



1 d/i 

e dx 
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The most relevant processes to determine the relaxation 
time T in electric conduction at low temperature are 
electron- impurity scatterings. Employing Born approx- 
imation, we can evaluate r by the chemical potential /i 
in the low temperature limit: 



(4) 



where njmp — Ni^^-p/V is the impurity density, m is the 
electron mass, h is Planck constant and uq is connected 
to the scattering potential V{r) as V{r) = uo J2i '^('^ ~ 
Ri) with the position of the random impurity Ri.^ 
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Fig. 1. The schematic picture of the electron transport: The spec- 
imen is directly connected with two equilibrium reservoirs. The 
reservoirs have chemical potential /x^ and /i^ respectively. The 
length of the specimen is L. 



From eqs. (3) and (4), the current density is given by 

dp 



3 
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where factor 2 is the spin degeneracy. Furthermore, we 
have to take into account the charge imbalances that ap- 
pear in the specimen. The electric field arises due to the 
local deviation of the electron density from its equilib- 
rium value and satisfies the Maxwell equation 



dE 
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where e is the permittivity and the unperturved electron 
density p{x) is 

dp 



p{x) 
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Now, we nondimensionalize the physical quantities 
(atomic units) as j <— jff {AnenY /m^e^ , ''^x ^ 
VxHAneo/e'^, x ^ xme^ /h^Aireo, p ^ /ift^(47reo)^/me'', 
E ^ Eh'^iATTeof/m'^e', n,^p ^ ni^p/i^ (47reo)V"^^e^ 
uq <~ ut^rri^e^ /K^A'neQ. Applying Sommerfeld expan- 
sion^^ to eq. (5), the dimensionless steady current density 
is given by 

2 1 , du, 
3^—-^—KE + i^), (8) 



where we notice that j is the negative value. Similarly, 
from eqs. (6) and (7), dimensionless Maxwell equation 
becomes 

dE 8 1 ^, dri^ 

dx \/6 ^r^o'^imp 



(9) 



with the relative permittivity er = e/eo- 

Let us consider the chemical potential inside the spec- 
imen. From eqs. (8) and (9), we obtain the ordinary dif- 
ferential equation of the chmeical potential: 



dx'^ jj? dx 
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Chemical potential profiles for the coefficient (a, b) 
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the chemical potential profile for the first order perturbative solu- 
tion with the corresponding coefficients. The ratio of the chemical 
potential ^il/ l^R is 1.1. The length of the specimen is L=2000. 
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Fig. 3. Top(A): Electrostatic potential profile for the coefficient 
(a,b)=(-10-*, -2 X 10-"') and (0,0). Bottom(B): Electrostatic 



potential profile for the coefficient (—10 



-10" 



). Line 1st are 



the electrostatic potential profile for the first order perturbative 
solution in eq. (13). All the remaining parameters are the same 
as in Fig. 2. 
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where dimensionless coefficients a and h are given by 



r= 0.0002 



2 '^o'i^impj 
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We now solve eq. (10) subject to the boundary conditions 
at /x(0) = iiL and /x(L) = fin. If we know the details of 
the material properties, the amount of the current den- 
sity under the chemical potential difference can be deter- 
mined. However, since the impurity potential Uq is hard 
to be measured in the experiment, we do not know the 
amount of the current density precisely. In this letter, we 
determine the spatial dependence of the chemical poten- 
tial as a function of a and b. 

Let us set the ratio iil/ ^J'R=^■^ with its normal- 
ized chemical potential ^l=0.0407, /Xi?=0.0370 and 
the normalized length L=2000 (/ZL=l.leV, fiR=1.0eV, 
L=106nm). Figure. 2 shows the spatial dependence of 
the chemical potental for the different coefficients a and 
b. In the finite current density, the chemical potential pro- 
files exhibit the strong nonlinearity. When the coefficient 
(a, b) is (— 10~*, —2 x 10~^°), the spatial variation of the 
chemical potential has a concave profile, while the chemi- 
cal potential profile with the coefficient (— 10~^, —10^^") 
shows a convex one. There exists the solutions of eq. (10) 
only when the value of b is smaller than 10~^ for any co- 
efficients a. 

The spatial dependence of the electrostatic potential 
in the specimen is of interest. The electrostatic potential 
is related to the electric field: (p = — J Edx. Using the 
electric field E = aiJ,~^ — d^/dx from eq. (8), the spa- 
tial dependence of the electrostatic potential can be de- 
termined. In Fig. 3(A), the electrostatic potential profile 
with the coefficient (— 10~*, —2 x 10~^°) shows a concave 
one. On the other hand, in Fig. 3(B), the electrostatic 
potential profile with the coefficient (—10^^, —10^^°) al- 
most increases linearly and shows a larger spatial varia- 
tion compared with the corresponding chemical potential 
profile in Fig. 2. 

If we treat the system with the electric filed, the elec- 
trostatic potential must be added to the chemical poten- 
tial for the potential of the electron. Now, we consider 
the dimensionless electrochemical potential /U = /x — in 
the specimen. We have shown that the chemical potential 
and the electrostatic potential have a strong nonlinear- 
ity for the finite current density, but the electrochemical 
potential profiles decrease almost linearly for any a and 
b in Fig. 4(A) and 4(B). If the gradient of the electro- 
chemical potential causes the electric current, it is rea- 
sonable that no electric current with the coefficient (0,0) 
shows the constant electrochemical potential profile in 
Fig. 4(A). Under the small external field, the linearity 
of the electrochemical potential profile corresponds the 
explanation of the electron transport by Ohmic law. 

These spatial dependence are reproducible by the per- 
turbation expansion. Now, we expand p., (j) and ft by the 
coefficient a: p, = p,o + ap,i + ■ ■ ■ , cf) = (f)Q + aipi + • • • and 
M = A*o+tiA'i + ' • ■ • Substituting these /x, </> and /x into eqs. 
(8) and (9), we obtain the zeroth order of a as /So = ci, 
(po = C2X + C3 and /io = /Jo + 4'o where the value of the 
each coefficient is given by ci = 0, C2 = — 1.85 x 10~^ 
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Fig. 4. Top(A): Electrochemical potential profile for the coeffi- 
cient (a,b)=(-10-*,-2 X 10-10) and (0,0). Bottom(B): Elec- 
trochemical potential profile for the coefficient (—10—^, —10~^^). 
Line 1st are the electrochemical potential profile for the first or- 
der perturbativc solution in eq. (12). All the remaining parame- 
ters are the same as in Fig. 2. 



and C3 = 0.0407 to recover the profile at (a, 6) = (0,0). 
Expanding until the first order of a, we obtain the elec- 
trochemical potential as 



— log(l -I- —x) 

C2 C3 ' 



(12) 



Similarly, the electrostatic potential is given by 

4 b 



= C2X + C3 ■ 



;{c2X + C3f^'^+diX + d2 



(13) 



where the coefficients di and ^2 are given by 

^1 = 7 \-Aiic2L + csf/' - - - log(l + ^L) 
L [3 aC2 C2 C3 

and 
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The first order perturbativc solution of the chemical po- 
tential is given by eqs. (12) and (13). For small electric 
currents, each perturbative solution plotted as the line 
1st almost correspond to the chemical potential profile 



in Fig. 2, the electrostatic potential profile in Fig. 3 and 
the electrochemical potential profile in Fig. 4 

Now, let us estimate the gap of the chemical poten- 
tial based on the approach by Nishino and Hayakawa.^ 
We first consider the distribution function at the connec- 
tion point. On the left reservoir side of the connection 
point, electrons with Px > obeying the equilibrium dis- 
tribution function /eg go through the connection point 
ballistically. In the same way, on the specimen side of 
the connection point, electrons with < obeying the 
nonequilibrium distribution function /„e go through the 
connection point ballistically. Therefore the dimension- 
less distribution function at the connection point is ap- 
proximately given by 



fcp{p) = 



/eq(P,Meq) Px > 



(14) 



fneip, /ine) Px < 0, 

where the equilibrium distribution function is 
f = I 

and the nonequilibrium distribution function as the first 
order of the current density is 

37r2 



(15) 



fne — fne 



fne{^-fne)Pv.f^n!^'j, (16) 



with /°g = l/(e''(^~''"«^-|-l). Theexpansionofthechem- 
ical potential in terms of the dimensionless current den- 
sity becomes 

l^ne = l^eq{i + dj), (17) 

where d is the gap parameter to be determined. Using 
eq. (14), we define the dimensionless current density at 
the connection point given by 

dp 



cp' 



(18) 



Since the electric current exists continuously at the con- 
nection point, we assume that J is equal to j in eq. 
(5). From these equations, the gap parameter is given by 
d = V37r72/v4,. Now, we consider the maximum amoimt 
of the gap at the left connection point. We have the so- 
lution of cq. (10), only when the coefficient b in eq. (11) 
is less than 10^®. If the specimen has the unit relative 
permittivity, the maximum amount of the gap is given 
by dj = -3.3 X 10"". Therefore, we confirm that the 
gap of the chemical potential at the connection point is 
negligible. 

Now, let us dicuss our result. Although we are inter- 
ested in the quantum effects of the electron transport, 
the essence of electrons is believed to be described by 
semi-classical treatments as we have discussed. Quantum 



effect may not be important for simple electron trans- 
ports, while such effects are essentially important in lo- 
calization problems. We have treated the system at the 
isothermal low temperature, but, it is straightforward to 
extend our analysis to the transport at high temperature 
cases, where the heat conduction is as important as the 
electron transport. In such a system, the gap of the ther- 
modynamical quantities at the connection point needs to 
be considered. The discussion of this effect is the subject 
of a subsequence paper. 

In summary, we have examined the electron transport 
processes which is driven by the chemical potential dif- 
ference based on BGK equation coupled with Maxwell 
equation at the isothermal low temperature. We deter- 
mine the spatial dependence of the chemical, electrostatic 
and electrochemical potential as functions of the current 
and the impurity potential. While the spatial dependence 
of the chemical potential and the electrostatic potential 
show the strong nonlinearity, the electrochemical poten- 
tial decreases almost linearly for any coefficient. The be- 
haviors of these thermodynamical quantities can be un- 
derstood by the simple perturbation theory. We also con- 
firm that the gap of the chemical potential at the con- 
nection point between the reservoir and the specimen is 
negligible. 
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